Symmetry energy from nuclear multifragmentation 
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The ratio of symmetry energy coefficient to temperature C sym /T is extracted from different pre- 
scriptions using the isotopic as well as the isobaric yield distributions obtained in different projectile 
fragmentation reactions. It is found that the values extracted from our theoretical calculation agree 
with those extracted from the experimental data but they differ very much from the input value 
of the symmetry energy used. The best possible way to deduce the value of the symmetry energy 
coefficient is to use the fragment yield at the breakup stage of the reaction and it is better to use 
the grand canonical model for the fragmentation analysis. This is because the formulas that are 
used for the deduction of the symmetry energy coefficient are all derived in the framework of the 
grand canonical ensemble which is valid only at the break-up (equilibrium) condition. The yield of 
"cold" fragments either from the theoretical models or from experiments when used for extraction 
of the symmetry energy coefficient using these prescriptions might lead to the wrong conclusion. 

PACS numbers: 25.70Mn, 25.70Pq 



Introduction:- Isospin-dependent phenomena in nu- 
clear physics has been an active area of research p| Q 
in recent years with the aim of enriching our knowl- 
edge about the symmetry term of the nuclear equation 
of state. The study of this quantity at different regimes 
of density and temperature is a hot topic [3-5] in the nu- 
clear physics community. This term plays an important 
role in areas of astrophysical interest such as the study 
of supernova explosions and the properties of neutron 
stars @. This also has significant influence in deciding 
the structure of neutron-rich and neutron-deficient nu- 
clei. The study of nuclear multifragmentation in heavy- 
ion reactions is an important tool for extracting infor- 
mation about the symmetry energy term, and this has 
created much interest in the nuclear physics community 
in recent years [7j-fl7| . 

In nuclear multifragmentation reactions, the neutron- 
proton composition of the break-up fragments is dictated 
by the asymmetry term of the equation of state and hence 
the study of the multifragmentation process allows one 
to obtain information about the symmetry term. Sta- 
tistical models (l8l - |20j that are simple and economic are 
very successful in predicting the outcome of nuclear mul- 
tifragmentation reactions. Different formulas [ij], l2lTj23T | 
have been proposed in the literature that connect the 
measurable fragment isotopic and isobaric observables of 
multifragmentation reactions to the symmetry energy of 
excited nuclei and these have been applied to the anal- 
ysis of heavy-ion collision data. These formulas have all 
been deduced using the grand canonical version of the nu- 
clear multifragmentation model assuming an equilibrium 
scenario for the break-up stage of the disintegrating sys- 
tem. They have been used to analyze experimental data 
from different projectile fragmentation as well as central 
collision reactions and the extracted values for the sym- 
metry e nerg y coefficient C sym ranges between 15 and 30 
MeV [7Hll| . It has been observed from model calcula- 
tions that the fragment distributions of the secondary 
fragments after evaporation vary from those at the pri- 



mary (breakup) stage and hence extraction of C sym val- 
ues from experimental yields that correspond to the 'cold' 
secondary fragments might lead to the wrong conclusion. 
This is primarily because the formulas that are used to 
extract the values of C sym from the experimental yields 
are all derived using the prescription of the grand canon- 
ical ensemble, which is valid only at the break-up (pri- 
mary stage) of the multifragmenting system. The yield of 
the fragments at the breakup stage of the reaction (before 
de-excitation by secondary decay) should only be used to 
extract information about C sym using the existing pre- 
scriptions. This is the main message of this paper and 
it has been established using projectile fragmentation re- 
actions. It has also been pointed out that theoretical 
simulations done by the grand canonical model arc only 
expected to reproduce the value of C sym used as input 
to the model. Calculations by canonical or any other 
ensemble that are otherwise better suited for describing 
intermediate energy nuclear reactions might lead to val- 
ues widely different from the input value of C sym used. 
Results from canonical and grand canonical ensemble 
[24| differ in general for finite nuclei [25[ . This work also 
presents a comparative analysis of the predictive power of 
the different existing formulas both at the primary stage 
and also after evaporation and their relative agreement 
with experimental data. 

Model description:- For the study of the nuclear sym- 
metry energy coefficient from projectile fragmentation 
reactions we use a model [2a - l28| that consists of three 
stages: (i) abrasion, (ii) multifragmentation, and (hi) 
evaporation. In heavy-ion collisions, if the beam energy 
is high enough then at a particular impact parameter 
three different regions are assumed to be formed: (i) the 
participant, (ii) the projectile like fragment (PLF) or pro- 
jectile spectator, and (hi) the target like fragment (TLF) 
or target spectator. Here we are interested in the frag- 
mentation of the PLF. 

In the abrasion part (first stage), we have calculated 
the PLF volume V s (b) at different impact parameters (b) 
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using straightline geometry. By knowing V s (b), we can 
determine the average number of neutrons and protons 
present in the projectile spectator (2(| by assuming the 
same neutron-to-proton ratio of projectile spectator and 
projectile. The impact parameter dependence of tem- 
perature is considered as T{b) = 7.5 — Ab[{A s {b) / Aq)] 
where A s (b) and Aq are mass numbers of projectile spec- 
tator and original projectile respectively. This tempera- 
ture profile successfully reproduces [13, HU the differen- 
tial charge distribution, total charge and mass distribu- 
tion, and isotopic distribution obtained in different ex- 
periments [29l432| . 

In the statistical models [l8l - [20| of nuclear multifrag- 
mentation, it is assumed that the hot nuclear system ex- 
pands and then, due to density fluctuations, it breaks 
up (multifragments) into composites of different sizes. 
Thermodynamic equilibrium is assumed to be reached 
at this expanded freeze-out configuration, where the in- 
teraction between the composites become unimportant 
(except for the long-range Coulomb force). The disin- 
tegration of hot nuclear system into various composites 
can be calculated using different statistical ensembles. In 
our model, in general canonical ensemble is used as it is 
more appropriate than the grand canonical for finite nu- 
clei calculations. Results from these ensembles differ in 
general though they are found to converge under certain 
conditions [25j for finite nuclei. The multifragmentation 
calculation of each PLF (second stage) formed at differ- 
ent impact parameters is done separately by using the 
canonical thermodynamical model (CTM) [l8[ which is 
based on equilibrium statistical mechanics and involves 
the calculation of partition functions. In this model the 
average number of composites with N neutrons and Z 
protons formed from a PLF N s , Z s is given by 
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where Qn s ,z s is the total partition function, which can 
be easily calculated using the recursion relations 33] . 

The grand canonical ensemble can also be used instead 
of the canonical, and the av erag e number of composites 
in this ensemble is given by [24j , 
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The neutron and proton chemical potentials ji n and fj, p 
are determined by the baryon and charge conservation 
conditions which amounts to solving the equations N s = 
J2Ne^ nN +^» z iJ N , z and Z s = '£Ze^ nN +^ z u; NtZ . 

In both the models, the partition function cjjv,.z of a 
composite having N neutrons and Z protons is a product 
of two parts and is given by 

u N ,z = ^(27rmT) 3 / 2 A 3 / 2 x z N , z (int) (3) 

The first part is due to the translational motion and the 
second part zjsr lZ (int) is the intrinsic partition function of 
the composite. V is the volume available for translational 



motion, which is the difference between the freeze-out 
volume and normal nuclear volume. For projectile frag- 
mentation calculation the freeze-out volume is assumed 
to be three times the normal nuclear volume. 

We now list the properties of the compos- 
ites used in this work. The proton and the 
neutron are fundamental building blocks, thus 
Zifi(int) = zo.i(mt) = 2, where 2 takes care of the 
spin degeneracy. For deuteron, triton, 3 He and 4 He we 
use ZN,z(int) — (2sjy,z + 1) exp[— /3En,z(9T)] where 
ft = l/T 7 EN,z(gr) is the ground-state energy of the 
composite and {2s ^.z + 1) is the experimental spin 
degeneracy of the ground state. Excited states for these 
very-low-mass nuclei are not included. For mass number 
A > 5 we use the liquid-drop formula. For nuclei in 
isolation, this reads 
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The expression includes the volume energy [Wo = 15.8 
MeV], the temperature dependent surface energy [cr(T) = 
a {(T 2 - T 2 )/(T 2 + T 2 )} 5 / 4 with a = 18.0 MeV and 
T c = 18.0 MeV], the Coulomb energy [a* = 0.31a c with 
a r = 0.72 MeV and Wigner-Seitz correction factor 0.31 
[19| ] and the symmetry energy (C sym — 23.5 MeV). The 
term 214 ( eo = 16.0 MeV) represents contribution from 
excited states since the composites are at a non-zero tem- 
perature. Different equations have been deduced using 
the grand canonical model which relates this symmetry 
energy coefficient C syrn to temperature as well as the 
source and fragment compositions. This is dealt with in 
the next section. 

The excited fragments produced after multifragmen- 
tation decay to their stable ground states. They can 
7-decay to shed energy but may also decay by light 
particle emission to lower mass nuclei. We include emis- 
sions of n,p, d, t, 3 He and 4 He. Particle-decay widths are 
obtained using Weisskopf's evaporation theory. Fission 
is also included as a de-excitation channel though for 
the nuclei of ^4<100 its role will be quite insignificant. 
The details of the evaporation stage (last stage) are 
described in Ref. [34l ]. 

Symmetry energy from different formulas:- There exist 
different formulas in the literature from which the sym- 
metry energy coefficient has been extracted. Here we 
make a short review of these existing formulas. 

The formula that connects the symmetry energy co- 
efficient with the isospin asymmetry of the source was 
first proposed in the framework of the Expanding evap- 
orating source (EES) model [211 ]. This is based on the 
grand canonical ensemble, and assuming thermodynamic 
equilibrium at the time of fragmentation of two systems 
having same charge Zq but different masses Aq±, Aq2 
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(Aqi < A02) at the same temperature T, C sym is given 
by 



a{Z)T 



(5) 



where a(Z) is the isoscaling parameter [2l], [35|, [36| of 
fragments having charge Z which can be measured from 
the ratio of the isotopic yields [37| ■ The suffix indicates 
the fragmenting source. This formula, which is referred 
to as the isoscaling (source) formula, has been extensively 
used on experimental data 0-fHl- 

The isoscaling (fragment) formula also derived [22| 
from the grand canonical ensemble assumes equilibrium 
at the breakup stage of two fragmenting sources in iden- 
tical thermodynamical states that differ in their isospin 
content (different masses Aq\, Aq2 but same charge Zq). 
This is given by 
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In this expression it is assumed that the isotopic distri- 
butions are essentially Gaussian and that the free ener- 
gies contain only bulk terms. Here, a(Z) is the isoscaling 
parameter of fragments having Z protons, < Ai > is the 
average mass number of a fragment of charge Z produced 
by source i=l (less neutron rich), 2 (more neutron rich), 
and T is the temperature of the decaying systems. This 
formula, which is similar in structure to the isoscaling 
(source) formula, connects the symmetry energy coeffi- 
cient to the isotopic composition of fragments instead of 
the isotopic composition of sources as in Eq. (5). 

An alternate expression (fluctuation formula) has been 
derived in Ref. [23l |38[ connecting the symmetry energy 
of a cluster of size A to the width of its isotopic distri- 
bution. Indeed, a Gaussian approximation on the grand- 
canonical expression for cluster yields gives 



Csym (^) 



AT 
2aj{Ay 



where aj (A) indicates the width of the isobaric distribu- 
tion of a cluster of size A and I = A — 2Z. 

According to the isobaric yield ratio method [hH - 
IT7L HH , obtained from the grand canonical expression for 
cluster yields, C sym can be expressed as 

AT , 



C sym (-^) 



-[lnR(3,l,A)-lnR(l,-l,A)} (8) 



where, R(I + 2, 1, A) is the yield ratio between two 
isobars differing by 2 units in isospin content as 
R(I + 2,I,A) = Y(I + 2,A)/Y(I,A). This equa- 
tion assumes the Coulomb terms in lnR(3, 1, A) and 
lnR(l,—l,A) are same and in both ratios mixing 
entropy terms are neglected [39j . 



Results:- The isotopic and isobaric yield distributions 
are related to the nuclear symmetry energy coefficient 
through different prescriptions, which have already been 
described briefly in the previous section. These relations 
are derived using the yields of the fragments obtained 
in the grand canonical ensemble. We first compare the 
symmetry energy coefficient obtained using the different 
prescriptions in both the canonical and the grand canon- 
ical ensemble. In this calculation we have used the sta- 
tistical models (canonical or grand canonical) in order 
to obtain the yields of the composites formed after mul- 
tifragmentation (second stage) of the hot single source. 
The sources used for the first two methods [Fig 1(a) and 
1(b)] are A 01 = 55 and A 02 = 60 and Z = 25, while for 
the later two methods [Fig 1(c) and 1(d)] where a sin- 
gle source is required the source Aq — 55 and Zq = 25 is 
used. The temperature used for the calculation is 5 MeV. 
These are the results for the primary fragments and no 
secondary decay is used for the de-excitation of the ex- 
cited primary fragments. There are some differences in 
the results from the canonical and the grand canonical 
ensemble, and these differences are almost same in all the 
methods used. In the canonical ensemble, the extracted 
symmetry energy coefficient changes with the fragment 
mass or charge and this variation is more or less same 
for all the four methods used for extraction of this coeffi- 
cient. The value of extracted symmetry energy coefficient 
varies between 25 and 50 MeV while the input symmetry 
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FIG. 1: (Color online) Symmetry energy coefficient obtained 
from canonical (red solid lines) and grand canonical model 
(black dotted lines) at T =5 MeV and input C sym =23.5 MeV. 
(a) and (b) represents the variation of extracted C sym with 
proton number Z by using Eq. (5) and Eq. (6) respectively 
from two sources having same Zo = 25 but different Ao = 
60 and 55 where as (c) and (d) indicates the variation of 
extracted Csym with mass number A calculated by using Eqs. 
(7) and (8) for Z = 25 and A = 55. 
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energy coefficient used is fixed at 23.5 MeV. 

In contrast, in the results from tfie grand canonical 
ensemble, C sym is independent of the fragment mass or 
charge. The value of C sym extracted from the grand 
canonical ensemble lies between 23 and 24 MeV and 
hence matches almost exactly with the input value used 
for the calculation. This difference in results between 
two ensembles is mainly because these formulas are all 
derived using the prescription of the grand canonical en- 
semble and hence when this ensemble is used to extract 
the value, the results agree with the input value. The 
results from the canonical ensemble deviate from that of 
the grand canonical ensemble in general for finite nuclei. 
Hence extraction of C sym leads to values that can differ 
widely from the input value used. 




FIG. 2: (Color online) Comparison of C a ym/T calculated 
by using (i) Eq. (5) (black dotted lines) (ii) Eq. (6) (red 
dashed lines) (iii) Eq. (7) (green solid lines) and (iv) Eq. (8) 
(blue dash dotted lines) for primary fragments in projectile 
fragmentation reactions Ni on Be (left panel) and Xe on Pb 
(right panel). For (i) and (ii), Csym/T is calculated from two 
sets of reactions 58 Ni on 9 Be and 64 Ni on 9 Be (left panel) and 
124 Xe on 208 Pb and 136 Xe on 208 Pb (right panel) and plotted 
against Z where as for (iii) and (iv), C sym /T is calculated for 
only neutron less reactions and plotted against A. 

In the next part, we show the results of the calculations 
from the projectile fragmentation model for two different 
reactions. Since for finite nuclei, the canonical ensemble 
is physically more acceptable than the grand canonical 
ensemble, in our model (as described previously) we have 
used the canonical ensemble for the fragmentation of the 
excited PLF. In this case we have used the reduced sym- 
metry energy coefficient which is nothing but the ratio 
of the symmetry energy coefficient to temperature i.e, 
Csym/T. Recent studies on this par ameter using differ- 



ent methods can be found in ref. [13Ml6| . This is mainly 
because it is difficult to estimate temperature both from 
experimental yields and from the model results. Hence 
it is better to express the results in the form of C sym /T 
instead of using only C sym since unambiguous extrac- 
tion of temperature is very difficult. In the isoscaling 
(source) formula [Eq. (5)], symmetry energy is related to 
the isoscaling parameter a and the Z/A value of the two 
sources. The other three formula depend solely on the 
properties of the fragments. In the isoscaling (fragment) 



formula [Eq. (6)], G sym /T depends on the isoscaling pa- 
rameter and on the Zj <A> values of the fragments. In 
the fluctuation formula [Eq. (7)], the coefficient depends 
on the width of the isobaric distribution of the fragments 
and their mass while in the isobaric yield ratio [Eq. (8)], 
it depends on the ratio of the fragment yields and their 
mass. In the break-up stage of the multifragmentation 
reaction, the yields of the primary fragments can be used 
to deduce the values of C sym /T from all the four for- 
mulas. The projectile fragmentation reactions involved 
are 58 Ni on 9 Be and 64 Ni on 9 Be at 140 MeV/n [HHH. 
From Eqs. (5) and (6), the variation of the reduced sym- 
metry energy coefficient C sym /T with fragment charge Z 
is shown in Fig. 2(a). In the same figure the results from 
Eq. (7) and (8) are shown as functions of the fragment 
mass A. The first two methods [Eq. (5) and (6)] depend 
on the isoscaling parameter which is calculated using two 
projectiles Ni 58 and Ni 64 having same value of Z. The 
fluctuation method [Eq. (7)] and the isobaric yield ratio 
method [Eq. (8)] depend on one source for the calcula- 
tion of the reduced symmetry energy coefficient and the 
projectile Ni 58 has been used. It is seen that the results 
from the primary fragments are close to each other for 
all four formula and C sym /T increases with the fragment 
mass or charge. Since these values are for the break-up 
stage before the final de-excitation, no comparison has 
been made to the experimental data. 

We have repeated the same calculation [Fig. 2(b)] for 
different projectiles (Xe 124 and Xe 136 ). The projectile 
fragmentation reactions involved were 124 Xe on 208 Pb 
and 136 Xe on 208 Pb at 1 GeV/n :32j. The beam energy 
is 1 GeV/n which is much higher than the previous value 
of 140 MeV/n [Hilj]. The trend of the results remains 
almost the same irrespective of the beam energy. The 
results from all four formulas are close to each other for 
the primary fragments (at the break-up stage) for a wide 
range of beam energies from 140 MeV/n to 1 GeV/n. 

In the next stage we investigated the results for the 
secondary fragments i.e, after the evaporation of the ex- 
cited primary fragments. In each figure, the results from 
the primary fragments are also shown for the sake of con- 
venience in comparison. In Fig. 3(a) we have shown the 
result for the isoscaling (source) formula [Eq. (5)]. Here 
Csym/T decreases after evaporation because the isoscal- 
ing parameter a decreases after evaporation at the tem- 
perature range used here while the denominator of the 
right-hand side of Eq. (5), which depends only on the 
source sizes remain unchanged. In both results from the 
model and from the experimental data, the reduced sym- 
metry energy coefficient does not depend very much on 
the fragment size. The numbers obtained from the exper- 
imental data are less than those obtained from our model 
for the projectile fragmentation. In Fig. 3(b) we have 
plotted the results for the isoscaling (fragment) formula 
[Eq. (6)]. In this case, C sym /T increases after secondary 
decay and the results more or less agree with those ex- 
tracted from the experimental data. The isoscaling 
(fragment) formula depends on the Z/<A> values of the 
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FIG. 3: (Color online) Variation of C sy m/T with atomic 
number Z calculated by using (a) Eq. (5) and (b) Eq. (6) for 
58 Ni on 9 Be and 64 Ni on 9 Be reactions. Panels (c) and (d) 
depict the variation of C sym jT with mass number A for 58 Ni 
on 9 Be reaction calculated from Eqs. (7) and (8) respectively. 
Experimental data (black squares) compared with theoretical 
results: primary fragments (blue dotted lines) and secondary 
fragments (red solid lines). 



fragments [see Eq. (6)]. For the less neutron rich source, 
this quantity does not change very much after evapora- 
tion but for the more neutron-rich nuclei, this quantity 
increases after evaporation since <A> decreases as the 
peak of the isotopic distribution shifts to the left (lower 
values of A) after evaporation. Hence the denominator of 
the right-hand side of Eq. (6) decreases. The isoscaling 
parameter a in the numerator also decreases after evap- 
oration but the denominator decreases much more and 
hence C sym /T increases after secondary decay as is seen 
from the Fig. 3(b). This is in contrast to the results from 
the isoscaling (source) formula where the denominator is 
independent of the property of the fragments. In Fig. 
3(c) we have plotted the results from the fluctuation for- 
mula [Eq. (7)]. Here also the trend of the results is same 
as in Fig. 3(b). C sym /T increases after secondary decay, 
and this is due to the fact that a 2 which is a measure 
of the width of the isobaric distribution, decreases after 
secondary decay and one can see from Eq. (7) that if this 
decreases then C sym /T will increase and this is exactly 
what happens as we see from Fig. 3(c). The experimental 
values obtained in this case also are quite close to those 
obtained from the model. In Fig. 3(d) we have plotted 
the results from the isobaric yield ratio method. In this 
case also the result is similar to that of the previous two 
cases and the reduced symmetry energy coefficient in- 
creases after evaporation. The numbers extracted from 
the experimental data are close to those from the theo- 
retical calculation. A similar trend of results from the 



isobaric yield ratio method in Antisymmetrized molecu- 
lar dynamics (AMD) model+GEMINI calculations is re- 
ported in [3]. In this case as seen from Eq. (8), C sym /T 
depends on the ratios of yields of different isobars. The 
reduced symmetry energy coefficient increases after evap- 
oration because the width of the isobaric distribution de- 
creases due to secondary decay, which results in decrease 
of the yield Y(-1,A) and Y(3,A) in Eq. (8) while the 
yield Y(1,A) remains almost unchanged. In the results 
from Eq. (6)-(8), C sym /T increases after evaporation, 
as compared to the results obtained from the primary 
fragments. 

We have repeated this calculation for another projec- 
tile fragmentation reaction (Xe on Pb) at 1 GeV/n [32| 
and this is shown in Fig. 4. The observations and re- 
sults are similar to those of the previous reaction in spite 
of the vast difference in the projectile energy and widely 
different target-projectile combination. 

All the four formulas are derived from the yields of 
the grand canonical ensemble and they hold good for 
the breakup stage of the reaction. Hence any attempt 
to deduce the value of the symmetry energy coefficient 
from the yields obtained after evaporation might lead to 
the wrong conclusion. Neither the experimental yields 
(which are the values after evaporation from the excited 
fragments) nor the yields of the secondary fragments 
from the model should be used to deduce the values of 
the symmetry energy coefficient. It might be possible 
to deduce the value of C sym from the break-up stage 
of the reaction, i.e., from the hot primary fragments, 
but since it is difficult to access this stage from an 
experimental point of view, no attempts are made to do 
such calculations. 
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FIG. 4: (Color online) Same as Fig.-3 except that here 
the projectile fragmentation reactions involved are 124 Xe and 
136 Xe on 208 Pb instead of 58 Ni and 64 Ni on 9 Be. 
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Summary and conclusion:- In this work we have com- 
pared the values of C sym obtained from the primary frag- 
ments for the canonical and the grand canonical ensem- 
bles for a single source at a fixed temperature using the 
four different prescriptions. The results from the canon- 
ical ensemble calculations show that the values obtained 
from the four formulas agree with each other. The results 
also display that for the canonical model, C sym varies 
with the fragment size and differs from the input C sym 
value which is equal to 23.5 MeV. In the grand canoni- 
cal model, this value is independent of the fragment size 
and is almost equal to the input value used. This is be- 
cause all the four formulas are deduced using the grand 
canonical ensemble at the break-up stage of multifrag- 
mentation. Hence it is possible to get back the value 
of the input symmetry energy coefficient using only the 
grand canonical model for calculating the yields of pri- 
mary fragments at the break-up stage. The results of 
canonical and grand canonical ensembles differ in gen- 
eral for finite nuclei 25J and hence in this case also the 
results from the two ensembles are different. 

We have also extracted the value of C sym /T from 
the yields of projectile fragmentation reactions using 
our model for projectile fragmentation. The results 
from the primary fragments are close to each other for 
all the four formulas used. In this model, where the 



canonical ensemble is used to calculate the yields of the 
hot primary fragments, the values of C sym /T obtained 
from the secondary fragments after evaporation are 
close to those obtained from the experimental yields 
but they differ from those obtained from the primary 
fragments and from the input value of C sym used in the 
model. The message we wish to convey is that to deduce 
the value of the symmetry energy coefficient using the 
prescriptions described earlier [Eqs. (5) to (8)], it is 
advisable to use the grand canonical model to obtain the 
yield of the fragments and one should use the yields at 
the break-up stage (primary fragments) of the reaction. 
The experimental yields (which are from the "cold" 
fragments) should not be used to deduce the value of 
the symmetry energy coefficient since the formulas used 
for the deduction are all valid at the equilibrium stage 
or the breakup stage of the reaction and secondary 
decay disturbs the equilibrium scenario of the breakup 
stage. Attempts to deduce the value of C sym from the 
secondary fragments or from the experimental yields 
might lead to values that are very different from the 
actual value. 

Acknowledgement:- The authors thank M. Mocko, 
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